#pragma once
#include "base/WrapperVtkmTypes.h"

namespace pde
{

class EquationSet
{
public:

	template<typename VariableType>
	VTKM_EXEC
	void FluxTerm(VariableType& uh)
	{
		VariableType conv_term[5];
		conv_term[0] = uh;
		conv_term[1] = uh;
		conv_term[2] = uh;
		conv_term[3] = uh;
		conv_term[4] = uh;
	}

public:
	EquationSet() {};
	~EquationSet() {};

private:
};

class InitalCondition
{
public:
	template<typename PointType, typename T>
	VTKM_EXEC
	void operator()(T time, const PointType & p) const
	{

	}

	InitalCondition() {};
	~InitalCondition() {};

private:

};


}